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Abstract 

In this paper, we test the consistency of Gamma Ray Bursts (GRBs) Data-set and Supernovae 

Union2 (SNU2) via the so-called multi- dimensional consistency test under the assumption that 

^~\ ' ACDM model is a potentially correct cosmological model. We find that the probes are inconsistent 

Q_c with 1.456o" and 85.47% in terms of probability. With this observation, it is concluded that GRBs 
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can be combined with SNU2 to constrain cosmological models. 
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I. INTRODUCTION 



The cosmic observations of type la supernovae (SN la) imply that our universe is under- 



going an accelerated expansion 



pe 1j 



2]. Furthermore, this implication was confirmed by the 



4] 



and large scale structure 



observations from cosmic microwave background radiation 3, 
5|, |6j. However, understanding the current accelerated expansion of our universe has become 
one of the most important issues of modern cosmology. In general, from the phenomeno- 
logical point of view, this late time accelerated expansion of our universe is due to possible 
modification of gravity theory at large scale or an exotic extra energy component, dubbed 
dark energy, which has negative pressure. 

To reveal the nature of the accelerated expansion or properties of dark energy, one needs 
more powerful cosmic probes. In the last decade, the data points of SN la (the current 
SNU2) have amounted to the number 557. However, the redshift range of SN la is relatively 
limited. Of course, higher redshift probes are useful to describe the evolution of our universe 
and to reveal the nature of late time accelerated expansion of our universe and properties 
of dark energy. The redshift of GRBs can extend to z ~ 8.1 or higher which makes it as a 
complementary cosmic probe to SN la. But before using GRBs to constrain cosmological 
models, the GRBs correlations, which relate cosmological models and intrinsic properties 
of GRBs, should be calibrated first. In general, the GRBs correlations can be written in a 
common form of y = a + bx where a and b are the calibrated parameters, x and y are re- 
ated to the intrinsic properties of GRBs and cosmological models, for the details please see 
7|. However, if one calibrates the GRBs correlations via a defined cosmological model, say 
ACDM model with Q m o = 0.27, the resulting distance moduli of GRBs are not independent 
of the input cosmology model. As a result, the obtained distance moduli can not be used to 
constrain any other cosmological models. The so-called circular problem will be committed 
once the above mentioned results are used to constrain any other cosmological model. Based 
on this point, the distance modulus derived by Schaefer |7j] can not be used to constrain any 
other cosmological models. So, new methods would be introduced to overcome this problem. 
Li, et. al [81] put the GRBs correlation and its cosmological model constraint together as a 
whole to fix the calibration parameters and to obtain the best fit values of the cosmolog- 
ical parameters in different cosmological models via Markov Chain Monte Carlo (MCMC) 
method. The lack of the GRBs calibration makes the GRBs weak to constrain cosmological 



models. In fact, the test of the correlations are needed to guarantee the consistency. In an 
alternative way, cosmography method was considered in 9j by parameterizing the luminosity 
distance di, in terms of deceleration parameter q Q , jerk j and snap Sq parameters. Clearly, 
the so-called circular problem is removed. However, this Taylor series method is limited 
when it is combined with higher redshift data point to constrain cosmological models. Liang 



et. al. 



10| used the low redshift SN la to calibrate the GRBs correlations and assumed 



the correlations were respected at high redshifts. Recently, this method was reconsidered 



by Wei ll|, Il2j . But, there would be some problems when GRBs are combined with other 



external data sets to constrain cosmological models. Wang presented a model-independent 



distance measurement f p (zi) (Eq. (TTUj) of this paper) from GRBs calibrated internally [131 ]. 
where zi are the redshifts of GRBs. The main points of Wang's method are that the resulted 
distance measurement f p (zi) is cosmological model in-dependent. The values of correlation 
parameters a and b are not used directly but the statistical errors of correlation parameters 
a a , <ib and systematic error <j sys are. This is because the la error bars of a, b and systematic 
error are almost the same in ACDM with different values of Q m o, though the values of a 
and b are different. Then, the cosmic constraint from GRBs is obtained in terms of a set of 
model-independent distance measurements. The merits of this method are listed as follows: 
(i) the constraint from GRBs is in a cosmological model independent way. So, it can be 
used to constrain other cosmological models, (ii) It is not calibrated by any other external 
data sets. It will not suffer any consistent problem when it is combined with other data sets 
as cosmic constraints, (iii) The cosmological model independent calibration is done first, 
(iv) Though the absolute calibration of GRBs is not known, the slopes of GRBs correlations 
can be used as cosmological constraints. Clearly, the drawback is clear that the constrained 
result is not tighter than the one calibrated by using SN la. But, if we have enough data 
points of GRBs, this problem will be overcome. Because the slopes of GRBs correlations 
are considered alone, this may make the GRBs not very powerful. 

Following the method proposed by Wang [13], Xu obtained N = 5 model-independent 
distances data sets and their covariance matrix by using 109 GRBs via Amati's E Pji — E iso 



correlation 23]. These five model- independent distances data points have been used to 
constrain cosmological model |24j . However, the consistency of the obtained five data sets 
via Amati's correlation with other cosmic probes must be checked to guarantee the reliability 
of GRBs. With this motivation, we will test the consistency or inconsistency of SNU2 with 



these five data sets derived from GRBs via the so-called multi- dimensional consistency test 
which will be reviewed briefly in section |H~C| for the details please see [251 ] . 

This paper is structured as follows. In section [XTJ, the SNU2, the five data sets derived 
from GRBs and the method to constrain dark energy model are presented. Also, the multi- 
dimensional consistency test is reviewed briefly. Section IIHI is the concluding remark. 



II. DATA-SETS AND METHOD 



A. Type la Supernovae 



Recently, SCP (Supernova Cosmo 
dataset which consists of 557 SN la [20] 



ogy Project) collaboration released their Union2 
The distance modulus fj,(z) is defined as 



H th (z) = 51og 10 [d L (z)] +/!o, 



(1) 



where d L {z) is the Hubble-free luminosity distance H dL(z)/c = H dA(z)(l + z) 2 /c, with 
H the Hubble constant, and no = 42.38 — 5 log 10 h through the re-normalized quantity h as 
Ho = lOOh km s~ Mpc~ . Where di,{z) is defined as 

dz' 



d L (z) = (1 + z)r(z), r(z) 



:Sinn 



VK\ 



E(z' 



(2) 



where E 2 (z) = H 2 (z)/Hq. Additionally, the observed distance moduli ii bs(zi) of SN la at 

Zi are 

^obs(zi) = m obs (zi) - M, (3) 

where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the parameters p s can be determined by a 
likelihood analysis, based on the calculation of 



x 2 ( Ps ,M>) = Yl 



{tlobs(Zi) - HthiPs, Zi)Y 



SN 



E 

SN 



(J 



{5 \og 10 [d L (p s , Zi)} - m obs (zi) + M'Y 



(4) 



where p s = {fi m o} denotes the model parameter and M' = /i + M is a nuisance parameter 
which includes the absolute magnitude and the parameter h. The nuisance parameter M' 



can be marginalized over analytically 2l| as 



X 2 (Ps 



-2 In / exp 

-oo 



-i x 2 ( Ps ,M') 



dM' 



resulting to 



? =A -* +ln (C- li 



(5) 



with 



SN 



A = S ^2 ( 51o Sio[^(Ps) z *)\ - m obs{zi)} ■ Covy 1 • {5\og w [d L (p s , Zj)] - m obs (zj)} , 

SN 

B = Yl Cov ^ ' i 5 lo Sio[^i(Pa, Zj)] - m obs (Zj)} , 

i 

SN 

C = ^Cov i : 1 , 



(6) 



where Cov.-,- is the inverse of covariance matrix with or without systematic errors. One can 

i fl 

find the details in Ref. [20J and the web site \j where the covariance matrix with or without 

systematic errors are included. Relation (Hj) has a minimum at the nuisance parameter value 

M' = B/C, which contains information of the values of h and M. Therefore, one can extract 

the values of h and M provided the knowledge of one of them. Finally, the expression 



X 2 SN ( Ps ,B/C) = A-(B 2 /C), 



(7) 



which coincides to Eq. (jHJ) apart from a constant, is often used in the likelihood analysis 



21 



22] . Thus in this case the results will not be affected by a flat M' distribution. It is 



worth noting that the results will be different with or without the systematic errors. In this 
work, all results are obtained with systematic errors. 



B. Gamma Ray Bursts 



Following [7], we consider the well-known Amati's E p j — Ei so correlation 14|, ll6Ml8| in 
GRBs, where E p j = E p ^ b s (l + z) is the cosmological rest-frame spectral peak energy, and 
E iso is the isotropic energy 

E lso = 47rd 2 L S bolo /(l + z) (8) 



1 http://supernova.lbl.gov/Union/ 



in which di and Sboio are the luminosity distance and the bolometric fruence of the GRBs 
respectively. Following j2|, we rewrite the Amati's relation as 



log 



Q, 



E 

-'-'ISO 

erg 



a + b log 



E„ 



300keV 



(9) 



In [131 ] . Wang defined a set of mo del- independent distance measurements {r p (zi)}: 



r p{z) 



r v\ z i) = ^T^' r v\ z ) 



[i + zy/*H 



r P (z ) 



-r(z) 



(10) 



where r(z) = dz,{z)/{\-\-z) is the comoving distance at redshift z, and z is the lowest GRBs 
redshift. Then, the cosmo logical model can be constrained by GRBs via 



Xgrbs(Ps) = [Ar p {zi)] ■ (Cov^ij ■ [Ar p ( Zi )}, 



Ar p (zi 



f d p ata (z t 



(11) 
(12) 



where r p (zi) is defined by Eq. (flOl) and (Ccw G ^ B )y,i, j = 1...N is the covariance matrix. In 
this way, the constraints from observational GRBs data are projected into the relative few 



quantities r p (z, 



p\^i)i 



1...N. 



Following the method proposed by Wang [13], Xu obtained N = 5 model-independent 
distances data sets and their covariance matrix by dividing 109 GRBs into five bins via Am- 
ati's E Pi i — E iso correlation [23j. The resulting model-independent distances and covariance 
matrix from 109 GRBs are shown below in Tab. [land Eq. (TT4|) . The {f p (zi)}(i = 1, ...,5) 





z 


fdata( z ) 


a(r p (z))+ 


a(r p (z)) 





0.0331 


1.0000 


- 


- 


1 


1.0000 


0.9320 


0.1711 


0.1720 


2 


2.0700 


0.9180 


0.1720 


0.1718 


3 


3.0000 


0.7795 


0.1630 


0.1629 


4 


4.0480 


0.7652 


0.1936 


0.1939 


5 


8.1000 


1.1475 


0.4297 


0.4389 



TABLE I. Distances measured from 109 GRBs via Amati's correlation with la upper and lower 
uncertainties [23]. zq = 0.0331 as lowest redshift was adopted. 



correlation matrix is given by 



(Cov GRB ) 



( 1.0000 0.7780 0.8095 0.6777 0.4661 \ 

0.7780 1.0000 0.7260 0.6712 0.3880 

0.8095 0.7260 1.0000 0.6046 0.5032 

0.6777 0.6712 0.6046 1.0000 0.1557 

\ 0.4661 0.3880 0.5032 0.1557 1.0000/ 



and the covariance matrix is given by 



(C0V GRB )ij = Cr{r p {z i ))(x{r p {z j )){CoVGRB)i3, 

where 

°"(*p(*t)) = a ( f p( Z i)Y 

a (r p (z,i)) = cr{fp{zi)) 
the a(f p (zi)) + and a(f p (zi))~ are the la errors listed in Tab. [H 



if f p (z)>f p (z) data ; 
, if f p (z) < f p (z) data , 



(13) 



(14) 

(15) 
(16) 



C. Method: Multi-dimensional Consistency Test 



In Ref. 25[, the multi-dimensional consistency test of probes was considered. If we 
consider M parameters and N probes, the method to test the consistency is to minimize 
the x 2 (^a) with respect to X a , 



M N 



i=l a,/3=l 



(17) 



where Xa is the best fit value returned from the ith probe with covariance matrix C^L and 
A a is a random point in cosmological space. In this case, the value of A a at the minimum 
X 2 is the best fit value. The goodness of fit is quantified by the value of x 2 i n the standard 
way, i.e., by checking the expectation value < Xmin >= v where v is the degrees of freedom. 
For example, in our case, we consider the ACDM model with Q m o as a free model parameter 
(M = 1) and two probes (N=2): SN and GRBs. So, the degrees of freedom [y = N — M) are 
2 — 1 = 1. The expectation value of Xmin would be 1. However, the value < Xmin >= v + B 
will be returned, where B > denotes the possible deviation from v. Then, the consistency 
can be concluded by the value of d a via the formula 



T(u/2,B/2)/T(u/2) = Erf(d a /V2). 



(18) 



The larger value of d a denotes better inconsistency between the probes. For example, a 
difference B = 9 tells us the two probes are inconsistent with 99.7% (3<r) in ACDM model. 
For convenience, we show the relation between B and probability in Fig. [TJ where v = 1 is 
adopted. 




FIG. 1. The relation between B and probability where v = 1 is adopted, where the red point 
denotes the probability 85.47% at B = 2.121. 

As described above, we firstly find the corresponding minimum Xmin values with SN and 
GRBs via Markov Chain Monte Carlo (MCMC) method. Our code is based on the publicly 
available CosmoMC package [26|. The results are shown in Tab. [Ill Via the formula 03 



we find the minimum value of x 2 (f2 m o) is 3.121 with the best fit value of Q m o = 0.287, where 
the covariance matrix C^l is given 



^ = *«(*%(*%* 



(19) 



here 



<W = <(AT, if A Q >A 
<(A*)=<(AT, if A Q <A 



% , 

Q ■ 



(20) 
(21) 



and the correlation between SN and GRBs is zero. However, the expected value of x 2 (^mo) 
should be v = 1. Then the returned deviation B is 2.121. From Eq. (TT8|) . one finds that the 
probes are inconsistent with d a = 1.456a and 85.47% in terms of probability. 

8 



Datasets 



SNU2 
GRB 



Parameters x^ in /d.0.f 



1 530.722/556 

1 3.0354/4 



n 



rnO 



n 97/ | +0.0386 
u - z '*-0.0358 

690 +0 ' 313 
U.DZU_ lg2 



TABLE II. The results of Xmin> ^mo with la regions are listed, where SN systematic errors are 
included, d.o.f denotes the degrees of freedom. 

III. CONCLUSION 

The consistency of SN la and GRBs is checked by the so-called multi- dimensional con- 
sistency test under the assumption that ACDM model is a correct cosmological model. We 
find that the inconsistency of SNU2 and GRBs is about 1.456er and 85.47% in terms of prob- 
ability. So, we can conclude that the five GRBs data points are consistent with SNU2 at 
the level above 1.456cr. These five GRBs data sets can be combined with SNU2 to constrain 
other cosmological models. 
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